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We study the CMB constraints on two Dark Energy models described by scalar fields with different 
Lagrangians, namely a Klein-Gordon and a Born-lnfeld field. The speed of sound of field fluctuations 
are different in these two theories, and therefore the predictions for CMB and structure formation are 
different. Employing the WMAP data on CMB, we make a likelihood analysis on a grid of theoretical 
models. We constrain the parameters of the models and compute the probability distribution 
functions for the equation of state. We show that the effect of the different sound speeds affects the 
low multipoles of CMB anisotropies, but is at most marginal for the class of models studied here. 
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Introduction 

Several observations now indicate 0, Q] that the uni- 
verse has been accelerating its expansion rate for the last 
5-10 Gy. Standard candidates for the "dark energy" 
(DE) which is responsible for this recent burst of ex- 
pansion are a cosmoiogicai constant and a scalar field 
IE BL IE IE which, for reasons unknown yet, started 
to dominate over the other types of matter just at our 
cosmoiogicai era. Unfortunately, DE seems to pose an 
even more formidable problem than that of dark matter 
to observational cosmology: whereas we still hope that 
the particles that make up dark matter will eventually 
be detected, and that incremental knowledge about the 
distribution of large-scale sctructure and galactic dynam- 
ics will eventually nail down the basic properties of dark 
matter and its relation to baryonic matter, no such hopes 
apply for DE, at least for now. Nevertheless, we must not 
shy away from trying to test models of DE, under penalty 
of theory running amok. 

One of the best opportunities to test DE is by looking 
at the anisotropies in the temperature and polarization 
of the cosmic microwave background radiation (CMB), 
which have been measured with exquisite accuracy by 
WMAP — and will be measured with even grander 
precision by the PLANCK mission [g. 

There are basically two ways through which DE can 
affect the CMB. The first is through its impact on the ex- 
pansion rate, which is determined by the equation of state 
of the DE component. The second is directly through 
the perturbations: DE, if it is not a plain cosmoiogicai 
constant, possesses small inhomogeneities which interact 
gravitationally with the inhomogeneities in baryons, dark 
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matter and relativistic matter. The physical properties of 
DE perturbations constitute additional ingredients which 
can impact the CMB anisotropies and LSS. 

Here we are interested in the role of the sound speed 
of DE perturbations. Whereas for the background evo- 
lution it is only necessary to specify the DE budget at 
the present time fixo, and its pressure/density ratio 
wx(t) = px{t)/px(t), in order to classify the perturba- 
tions sector we need to specify the pressure perturbations 
as well 0: 

Spx = c x Spx + 3H(1 + wx) 6 -^f (c\ - , (1) 

where px , Spx , c x and Ox are respectively the DE den- 
sity, density perturbation, sound speed and velocity po- 
tential. The case for a perfect fluid DE model in which 
Spx — c x 5px has been analyzed in 0] and compared 
with the WMAP data in In this paper we focus on 
the scalar field case in which the pressure perturbation 
is not simply related to the density perturbation. 

The standard quintessence scenario (a canonical scalar 
field described by a Klein-Gordon Lagrangian) is a min- 
imal modification to ACDM, and in that particular case 
c x = 1. All the information contained in the scalar po- 
tential V(tp) is encoded in wx in that case. 

The possibility of having c\ ^ 1 in the context of 
scalar field theories by considering a Lagrangian with 
a non-standard kinetic term is not new. If such non- 
canonical scalar field plays the role of an inflaton, the 
predictions for scalar and tensor perturbations are differ- 
ent with respect to the usual (canonical) scenario [l2j ■ In 
the context of DE, K-essence [j| was proposed (see also 
0), even if it has the unpleasant feature of having fluc- 
tuations which can travel with speeds faster than that of 
light. The imprints of K-essence on CMB anisotropies 
have already been studied elsewhere |E El , but not yet 
in a statistic way. Here we focus on a different scalar field 
theory and we perform a statistical study of the predic- 
tions 113 • 
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To be concrete we study a Klein-Gordon (KG hence- 
forth) Lagrangian density: 
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(2) 



and compare it to a Born-Infeld 15] (BI henceforth) La- 
grangian 
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(where we have assumed a signature ( — h H — h) for the 
metric and we have introduced the mass scale M). This 
Lagrangian can be thought as a field theory generaliza- 
tion of the Lagrangian of a relativistic particle [l6j . For 
d^cjjd^/M 4 « 1, by taking the first term of a Tay- 
lor expansion of Eq. (|3J), we get a Lagrangian which is 
equivalent to the KG one by redefining the field (except 
for problems with the invertibility of this redefinition). 
Only by taking the full Lagrangian ©, the two theo- 
ries are different. The BI was recently revisited in con- 
nection with string theory, since it seems to represent a 
low-energy effective theory of D-branes and open strings 
[ttL ITsf . In the case of constant potential, the model re- 
duces to a type of hydrodynamical matter known as the 
Chaplygin Gas 0] — whose predictions were studied in 

[UGl. While for the KG La 

grangian the sound speed of 
fluctuations is 1 independently of the potential, in the BI 
case c x — —wx > always, regardless of the potential. 

For both theories we will assume an inverse power-law 
potential: 
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where p > and / stands for ip or </> in the canonical 
or non-canonical cases, respectively. Positive powers of 
the field (p < 0) yield models which, if not finely tuned, 
do not lead to acceleration. In the following we split the 
fields as / = fo(t)+Sf(t, x), as usually done in cosmology. 



Equation of state versus sound speed 

In the canonical case, the scalar field model with po- 
tential 10} is known as the Ratra-Peebles model of dark 
energy 3]. In the BI case it was shown in that this 
model can lead to an accelerating regime when p < 2. 
With the potential given by Eq. , both theories behave 
very similarly at the level of the background. This is so 
because both scalar fields possess fluid-like attractor solu- 
tions when the background energy density is dominated 
by a perfect fluid such as dust (wf = 0) or radiation 
(wf = 1/3). For the Ratra-Peebles model the equation 
of state of the attractor during a fluid-dominated period 
is 201: 
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FIG. 1: Upper panel: background quantities for two similar 
Born-Infeld (BI, dashed lines) and Klein-Gordon (KG, solid 
lines) dark energy models as a function of redshift. The upper 
curves denote Q r and Q, m , and the lower, nearly superimposed 
curves denote the equations of state for the two models. The 
cosmological parameters are fit = 0.04, Q m = 0.26, h 
0.7'.'. and either p = 3/2 (BI) or p = 6 (KG.) Lower panel: 
temperature angular power spectra for KG model (solid line), 
BI model (dashed line) and a fiducial ACDM (dotted line), all 
normalized at the location of their first doppler peak. Only 
for comparison we also show the (wrong) power spectrum of 
the BI model obtained without considering the scalar field 
perturbations (dot-dashed line). 



where ifo (t) is the homogeneous value of the scalar field. 
For the BI model the equation of state during a fluid- 
dominated period is [2j: 
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In particular, for p <^ 1, w v ~ w^, while tracking is more 
and more accurate for p >> 1 in the KG case and for 
p ~ 2 in the BI case. 

The main difference between them lies in the behaviour 
of their perturbations. The canonical scalar field pertur- 
bations difi^x, t) obey the equation of motion: 
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where a is the scale factor and h is the synchronous gauge 
metric perturbation. The BI scalar perturbation, on the 
other hand, obeys the equation: 
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+ M i (\ogV)^8<i>=--hfo. (8) 

Eqs. of motion Q and JSJ determine how perturbations 
behave: while the canonical scalar field perturbations 
have a speed of sound c£ = 1, in the BI case the speed of 

sound is c| = — w$ — 1 — </>q/M 4 . Since c x determines 
how fast fluctuations dissipate, a lower sound speed in- 
creases the phase space of modes which are Jeans unsta- 
ble [3] [to make contact with Eq. JTJ, 9 X = k 2 5f/f, 
where 5f is the field perturbation.] Hence, the spectra 
of anisotropics predicted by the two models are slightly 
different because of the different Jeans scales in the DE 
sector - see Fig. 1. 

A key feature of these models is that not only their 
backgrounds are well described by a simple attractor: 
their perturbations are funnelled to an attractor as well 
plf. This means that the issue of initial conditions for 
the ba ckgr ound field and the perturbations is partially 
solved 28|. If this was not true, initial conditions would 
constitute additional free parameters in the models. 

CMB phenomenology and likelihood analysis 

In order to compute the CMB anisotropies we em- 
ployed a modified version of the code CMB FAST v. 4.1 
[22| . adding scalar fields to the system of equations. 

We take the following vector of free parameters: 

a = {H,Q, m ,Q, b ,n s ,m,p;A} , (9) 

where H denotes the present value of the Hubble constant 
in Km s _1 Mpc -1 , Q, rn is the amount of dark matter, f2f, 
the amount of baryons, n s is the scalar spectral index, m 
is the scalar field mass scale in the potential (@J and p is 
the power of the scalar potential. As always, there is a 
free parameter (A) related to the (unobservable) overall 
normalization of the CMB spectrum. Since we consider 
only flat models, the density in the dark energy compo- 
nent, fix, is fixed once the other cosmological parameters 
are determined. We have ommitted the BI scale M from 
the space of parameters since it can be absorbed by a 
redefinition of m. 

We have evaluated the likelihood function C(Ci\a) for 
a grid of roughly 370,000 models using the code described 
in [23|. In order to get the posterior probability distri- 
bution function (p.d.f) for a given parameter, we must 
marginalize the likelihood function over the remaining 
ones. Since C(Ci\a) is not a p.d.f in the usual sense, we 
use Bayes' Theorem: 
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FIG. 2: Region of [H, Q m ] probed by our simulations. 

< tt b < 0.069, 0.95 < n s < 1.09. In the matter era 
(u>p = 0), the power p of the potential can be substituted 
by ujf because of Eqs. ©-© , so we chose a grid that was 
equally spaced in u) t rather thanp, with —0.6 < u>f < — 1. 
Regarding f2 m and H , it was numerically more suitable to 
use a certain combination of the two parameters which 
effectively explored a region of parameter space whose 
shape is shown in Fig. 2. 

We can study the parameters H and u>f (or, equiv- 
alently, p) jointly by marginalizing against (i.e., inte- 
grating out) all other parameters according to (ftTHl . thus 
obtaining the posteriori marginalized p.d.f. P m (H,ujf). 
This joint p.d.f for H and ojj is presented in Fig. 3 for 
the models under scrutiny. The contours correspond to 
1, 2 and 3cr levels. For the Ratra-Peebles (KG) case, the 
joint p.d.f for H and ujf agrees with [24|. 

We are mainly interested in ujf, so we marginalize 
against H as well to obtain the posteriori p.d.f. for the 
equation of state of dark energy. The p.d.f's, shown in 
Fig. 4, have been smoothed through a Bezier spline so 
that the effects of having an imperfect grid are not over- 
estimated. 

The limits on the equation of state of the dark energy 
component are, for the KG model: 

lu v > -0.84 (la), -0.74 (2a), -0.68 (3a) , (11) 

where la, 2a and 3a correspond to 68.3% C.L., 95.4% 
and 99.7% C.L., respectively. For the BI model the limits 
on its equation of state are: 

w > -0.87(lcr), -0.76 (2a), -0.69 (3a) . (12) 



P(ai) oc J L(Ci\a)Y[V(a j )da j , (10) 

where V(ctj) is the prior p.d.f for aj. 

For the parameters in @ we set an equally spaced 
grid with the following (flat) priors: 64 < H < 80 , 0.036 



Conclusions 

We have tested two scalar field models of dark en- 
ergy against the WMAP data on CMB anisotropies. For 
simplicity we have taken the same potential for both 
Lagrangians. Our aim was to compare models which 



4 



-1 -0.9 -0.8 -0.7 



74 



72 




FIG. 3: Joint posteriori p.d.f. for H and u f from WMAP 
data, for Ratra-Peebles (upper panel) and Born-Infeld (lower 
panel) models. The shaded areas correspond to 1, 2 and 3<r 
levels. 

had nearly identical backgrounds, but whose perturba- 
tive sectors behave differently. We did this by compar- 
ing two scalar field models, a Klein-Gordon (the usual 
quintessence) and a Born-Infeld scalar, both of which 
have very similar attractor solutions during radiation- 
and matter-domination. 

From a fundamental perspective, we have shown that a 
BI scalar field can easily play the role of DE. As for being 
a CDM candidate 25] the non-linear stage must still be 
carefully analyzed |2(| . We stress that a very interesting 
feature of BI theories is that the BI scalar can act as 
dust and drive the universe into acceleration, where the 
trigger to the accelerated phase is the transition of the 
background from radiation-dominated to dust-dominated 
at z eq ~ 10 4 0. 

For the Ratra-Peebles model (a canonical scalar field 
model) we find an allowed range of parameters p < 0.38 
at la, corresponding to an equation of state during 
the attractor regime in the matter-dominated period of 
— 1 < Wm < —0.84. The Born-Infeld model, on the other 
hand, is more tightly constrained: we find p < 0.27 at 



lcr, corresponding to an equation of state of the attrac- 
tor regime of — 1 < < —0.87. For the Born-Infeld 
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FIG. 4: Probability distribution function as a function of 
the equation of state for both models: KG (solid line) and BI 
(dashed line.) 



model, this also implies a limit for the sound speed of 
dark energy > 0.87 at the la level. 

The effect of the perturbations was not dramatic in the 
BI case, since we have c x = —wx (compared to c 2 x = 1 
for quintessence). However, in the cases studied here, the 
speed of sound affects interestingly the low-£ multipolcs. 
In order to have a DE model which fits observations well, 
one has to consider w x — 0.5 and therefore the sound 
speed for BI fluctuations is not sufficiently different from 
quintessence. Because of cosmic variance and high error 
bars, low £'s have a diminished influence on the Likeli- 
hood, which means that the region where the different 
perturbative behaviors show up more conspicuously are 
not well represented in the Likelihood analysis. Never- 
theless, changing the speed of sound of dark energy by a 
factor of order 10% can have an impact (from the point 
of view of the Likelihoods) of the same order of magni- 
tude as changing the other cosmological parameters by 
a few percent, which is the precision target of future ex- 
periments such as PLANCK @. 

This reinforces the notion that only dark energy 
models which possess a highly unusual perturbative be- 
havior, such as certain models of k-essence models 
with c x ~ (and unrelated to wx) or dark energy per- 
fect fluids 0, ^3 (where the pressure perturbation in 
Eq. Q depends only on the density perturbation) can 
have a large impact on the CMB. 
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